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We study the temperature-dependent conductivity o~(T) 
and spin susceptibility x(^) °f the two-dimensional disor- 
dered Hubbard model. Calculations of the current-current 
correlation function using the Determinant Quantum Monte 
Carlo method show that repulsion between electrons can sig- 
nificantly enhance the conductivity, and at low temperatures 
change the sign of da/dT from positive (insulating behav- 
ior) to negative (conducting behavior). This result suggests 
the possibility of a metallic phase, and consequently a metal- 
insulator transition, in a two-dimensional microscopic model 
containing both repulsive interactions and disorder. The 
metallic phase is a non-Fermi liquid with local moments as 
deduced from a Curie-like temperature dependence of x{T)- 



When electrons arc confined to two spatial dimensions 
in a disordered environment, common understanding un- 
til recently was that the electronic states would always be 
localized and the system would therefore be an insulator. 
This idea is based on the scaling theory of localization 
for non-interacting electrons and corroborated by subse- 
quent studies using renormalization group (RG) meth- 
ods The scaling theory highlights the importance of 
the number of spatial dimensions and demonstrates that 
while in three dimensions for non-interacting electrons 
there exists a transition from a metal to an Anderson in- 
sulator upon increasing the amount of disorder, a similar 
metal-insulator transition (MIT) is not possible in two 
dimensions. 

The inclusion of interactions into the theory has been 
problematic, certainly when both disorder and interac- 
tions are strong and perturbative approaches break down. 
Following the scaling theory the effect of weak interac- 
tions in the presence of weak disorder was studied by dia- 
grammatic techniques and found to increase the tendency 
to localize H . Subsequent perturbative RG calculations, 
including both electron-electron interactions and disor- 
der, found indications of metallic behavior, but also, for 
the case without a magnetic field or magnetic impurities, 
found runaway flows to strong coupling outside the con- 
trolled perturbative regime and therefore were not con- 
clusive The results of such approaches therefore 



have not changed the widely held opinion that in two 
dimensions (2D) the MIT does not occur. 

The situation changed dramatically with the recent 
transport experiments on effectively 2D electron systems 
in silicon metal-oxide-semiconductor field-effect transis- 
tors (MOSFETs) which have provided surprising evi- 
dence that a MIT can indeed occur in 2D §. In these 
experiments the temperature dependence of the conduc- 
tivity <7dc changes from that typical of an insulator (de- 
crease of CTdc upon lowering T) at lower density to that 
typical of a conductor (increase of Ode upon lowering T) 
as the density is increased above a critical density. The 
fact that the data can be scaled onto two curves (one 
for the metal, one for the insulator) is seen as evidence 
for the occurrence of a quantum phase transition with 
carrier density n as the tuning parameter. The possibil- 
ity of such a transition has stimulated a large number of 
further experimental |?],|| and also theoretical investiga- 
tions pjl0|, including proposals that a superconducting 
state is involved |TT| . Explanations in terms of trapping 
of electrons at impurities, i.e. not requiring a quantum 
phase transition have also been put forward p2| . While 
there is no definitive explanation of the phenomena yet, 
it is likely that electron-electron interactions play an im- 
portant role. 

The central question motivated by the experiments is: 
Can electron-electron interactions enhance the conduc- 
tivity of a 2D disordered electron system, and possibly 
lead to a conducting phase and a metal-insulator tran- 
sition? It is this question that we address by studying 
the disordered Hubbard model which contains both rele- 
vant ingredients: interactions and disorder. While the 
Hubbard model does not include the long range nature 
of the Coulomb repulsion, studying the simpler model 
of screened interactions is an important first step in an- 
swering the central qualitative question posed above. We 
use Quantum Monte Carlo simulation techniques which 
enable us to avoid the limitations of perturbative ap- 
proaches (while of course being confronted with others). 
We mention that recent studies using very different tech- 
niques from ours have indicated that interactions may en- 
hance conductivity: two interacting particles instead of 
one in a random potential has a delocalizing effect [fl3|| , 
and weak Coulomb interactions were found to increase 
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the conductance of spinless electrons in (small) strongly 
disordered systems Q. 

The disordered Hubbard model that we study is de- 
fined by: 



ij c ia C jcr 



3,<? 



(1) 



where Cj a is the annihilation operator for an electron at 
site j with spin a. tij is the nearest neighbor hopping in- 
tegral, U is the on-site repulsion between electrons of op- 
posite spin, /i the chemical potential, and n 3a = 
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the occupation number operator. Disorder is introduced 
by taking the hopping parameters £y from a probability 
distribution P(*y) = 1/A for 4 y € [1 - A/2, 1 + A/2], 
and zero otherwise. A is a measure for the strength of 



the disorder | 
We use t 
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re Determinant Quantum Monte Carlo 
(QMC) method, which has been applied extensively to 
the Hubbard model without disorder |l6} ]. While dis- 
order and interaction can be varied in a controlled way 
and strong interaction is treatable, QMC is limited in 
the size of the lattice, and the sign problem restricts the 
temperatures which can be studied. The sign problem is 
minimized by choosing off-diagonal rather than diagonal 
disorder, as at least at half-filling ((n) = 1) there is no 
sign problem in the former case, and consequently simula- 
tions can be pushed to significantly lower temperatures. 
For results away from half filling we choose (n) = 0.5 
where the sign problem is less severe compared to other 
densities J16[ . Also, interestingly, the sign problem is 
reduced in the presence of disorder yl| . 

The quantity of immediate interest when studying a 
possible metal-insulator transition is the conductivity 
and especially its T-dependence. By the fluctuation- 
dissipation theorem o~dc is related to the zero-frequency 
limit of the current-current correlation function. A com- 
plication of the QMC simulations is that the correla- 
tion functions are obtained as a function of imaginary 
time. To avoid a numerical analytic continuation pro- 
cedure to obtain frequency-dependent quantities, which 
would require Monte Carlo data of higher accuracy than 
produced in the present study, we employ an approx- 
imation that was used and tested before in studies of 
the superconductor-insulator transition in the attractive 
Hubbard model This approximation is valid when 

the temperature is smaller than an appropriate energy 
scale in the problem. Additional checks and applicabil- 
ity to the present problem are discussed below. The ap- 
proximation allows (Tdc to be computed directly from the 
wavevector q- and imaginary time r-dependent current- 
current correlation function A xx (q, t) : 
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<Jdc = — A ra (q = 0,r = /3/2) 

7T 



(2) 



Here = 1/T, A xx (q,r) = (j x (q, r)j x (-q, 0)), and 
jx (q, t) the q, r-dependent current in the x-direction, is 



the Fourier transform of j x (£) = i U+x,t{c\ +x>cr c 

4a C i+x,a) ■ ( SeC als ° Ref - II)- 
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FIG. 1. Conductivity Ode as a function of temperature 
T for various values of disorder strength A at U = 4 for (a) 
half-filling ((n) = 1) and (b) (n) — 0.5. Calculations are 
performed on an 8 x 8 square lattice; data points are averages 
over 4 realizations for a given disorder strength. 

As a test for our conductivity formula (3) we first 
present results in Fig. 1(a) for (Jdc{T) at half-filling for 
U = 4 and various disorder strengths A. The behav- 
ior of the conductivity shows that as the temperature 
is lowered below a characteristic gap energy, the high T 
"metallic" behavior crosses over to the expected low T 
Mott insulating behavior for all A, thereby providing a 
reassuring check of formula (2) and our numerics. 

In Fig. 1(b), we show (7d c (T) for a range of disorder 
strengths at density (n) = 0.5 and U = 4. The figure 
displays a striking indication of a change from metallic 
behavior at low disorder to insulating behavior above a 
critical disorder strength, A c ~ 2.7. If this persists to 
T = and in the thermodynamic limit, it would describe 
a ground state metal-insulator transition driven by 



2 



disorder. 
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FIG. 2. Conductivity a^ c as a function of temperature 
T comparing U = 4 and [7 = at (n) = 0.5 and disorder 
strength A = 2.0. Data points are averages over many re- 
alizations for this disorder strength (see text). Error bars 
are determined by the disorder averaging and not the Monte 
Carlo simulation. 



system we note that at lower values of A, where the local- 
ization length exceeds the lattice size, <7o shows "metal- 
lic" behavior which is diminished upon turning on the 
interactions Based on our analysis above, we would 
predict that at low enough T and large enough lattice 
size, the conductivity curves for the non-interacting ero 
and interacting a cross with a > aX sufficiently low T. 
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In order to obtain a more precise understanding of the 
role of interactions on the conductivity, we compare in 
Fig. 2 the results of Fig. 1(b) with the disordered non- 
interacting (To fl9|| . The comparison is made at strong 
enough disorder A = 2.0 such that the localization length 
is less than the lattice size and the non-interacting system 
is therefore insulating with da$/dT > at low T. Inter- 
actions are found to have a profound effect on the conduc- 
tivity: in the high-temperature "metallic" region, inter- 
actions slightly reduce a compared to the non-interacting 
cto behavior. On the other hand in the low-temperature 
"insulating" region of oq the data shows that upon turn- 
ing on the Hubbard interaction the behavior is com- 
pletely changed with da/dT < 0, characteristic of metal- 
lic behavior. This is the regime of interest for the MIT. 

In order to ascertain that the phase produced by re- 
pulsive interactions at low T is not an insulating phase 
with a localization length larger than the system size but 
a true metallic phase we have studied the conductivity 
response for varying lattice sizes. We find a markedly 
different size dependence for the {7 = insulator and the 
{7 = 4 metal, resulting in a confirmation of the picture 
given above. For {7 = 0, the conductivity on a larger 
(12 x 12) system is lower than that on a smaller (8 x 8) 
system (see Fig. 2), consistent with insulating behavior 
in the thermodynamic limit, whereas for {7 = 4 the con- 
ductivity on the larger (8 x 8) system is higher than that 
on the smaller (4 x 4) system (data not shown), indica- 
tive of metallic behavior. Thus the enhancement of the 
conductivity by repulsive interactions becomes more pro- 
nounced with increased lattice size |2(J] . 

Concerning finite-size effects for the non-interacting 



FIG. 3. Spin susceptibility x as a function of temperature 
T at (n) = 0.5 comparing interaction strengths U = 0, 2, 4 
and disorder strengths A = 2, 4. Calculations are performed 
on 8 x 8 square lattices; error bars are from disorder averages 
over up to 8 realizations. 

To obtain information on the spin dynamics of the 
electrons and because it is a quantity often discussed in 
connection with the localization transition, we also com- 
pute the spin susceptibility ^asa function of T (through 
x(T) — PSo(T) where 5*o is the magnetic structure fac- 
tor at wavevector q = 0). Fig. 3 shows two things: 
1) x{T) is enhanced by interactions with respect to the 
non-interacting case (at fixed disorder strength), and 2) 
starts to diverge when T is lowered, both on the metallic 
(A = 2) and insulating (A = 4) sides of the alleged tran- 
sition. This is in agreement with experimental and the- 
oretical work on phosphorus-doped silicon, where a (3D) 
MIT is known to occur and the behavior is explained by 
the existence of local moments p2| , and also with dia- 
grammatic work on 2D disordered, interacting systems 
©■ 

In order to definitively establish the existence of a pos- 
sible quantum phase transition in the disordered Hub- 
bard model requires: (i) Extending the present data at 
T = 0.1 = W/80, where W is the non- interacting band- 
width, to lower T, which is however difficult because of 
the sign problem, (ii) A more detailed analysis of the 
scaling behavior in both linear dimension and some scaled 
temperature, (iii) A more accurate analytic continuation 
procedure to extract the conductivity. The condition for 
the validity of the approximate formula (|^) for o~d c (T), 
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requires that T be less than an appropriate energy scale 
which is fulfilled within the two phases, but breaks down 
close to a quantum phase transition where the energy 
scale vanishes. 

In summary, we have studied the temperature- 
dependent conductivity cr(T) and spin susceptibility 
x(T) of a model for two-dimensional electrons containing 
both disorder and interactions. We find that the Hub- 
bard repulsion can enhance the conductivity and lead to 
a clear change in sign of da/dT. More significantly, from 
a finite size scaling analysis we demonstrate that repul- 
sive interactions can drive the system from one phase to 
a different phase. We find that a(T) has the opposite 
behavior as a function of system size in the two phases 
indicating that the transition is from a localized insulat- 
ing to an extended metallic phase. The x(T) data further 
suggests the formation of an unusual metal, a non- Fermi 
liquid with local moments. While the simplicity of the 
model we study prevents any quantitative connection to 
recent experiments on Si-MOSFETs, there is neverthe- 
less an interesting qualitative similarity between Fig. 1(b) 
and the experiments. Varying the disorder strength A 
at fixed carrier density (n), as in our calculations, can 
be thought of as equivalent to varying carrier density 
at fixed disorder strength, as in experiments, since in 
a metal-insulator transition one expects no qualitative 
difference between tuning the mobility edge through the 
Fermi energy (by varying A) and vice-versa (by varying 
(n)). Our work then suggests that electron-electron in- 
teraction induced conductivity plays a key role in the 2D 
metal-insulator transition. 
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